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ABSTRACT 

We present an analysis of two strongly lensed galaxies at z = 6.56 and z ^ 7 for which 
multi-band photometric and spectroscopic observations are available. For one source 
the data include recent HST and Spitzer observations. Using an SED fitting technique 
considering a number of parameters (various libraries of empirical and theoretical 
^ . template spectra, variable extinction and extinction laws) we attempt to constrain 

' the properties of their stellar populations (age, star formation history, mass) and their 

OO . intrinsic Lya emission. The following main results are obtained for the individual 

' galaxies: 
\0 ■ 

■ • Triple arc in Abell 2218, probable z ~ 7 galaixy: The most likely redshift 
ly^ \ of this source is z ~ 6.0-7.2 taking into account both our photometric determination 

■ and lensing considerations. SED fits indicate generally a low extinction {E{B — V) < 

0.05) but do not strongly constrain the star formation (SF) history. Best fits have 
typical ages of ~ 3 to 400 Myr. A reasonable maximum age of (250-650) Myr 
(1 a interval) can be estimated . However, the apparent 4000 A break observed 

O \ recently from combination of IRAC/Spitzer and HST observations, can also equally 

well be reproduced with the template of a young (~ 3-5 Myr) burst where strong 
^ ' restframe optical emission lines enhance the 3.6 and 4.5 /im fluxes. The estimated 

^ ' SFR is typically ~ 1 M© yr"i for a Salpeter IMF from 1-100 Mq, in agreement 

with previous estimates. The unknowns on the age and star formation history could 
easily explain the apparent absence of Lya in this galaxy. 
, • Abell 370 HCM 6A: The available Lya and continuum observations indicate 

' basically two possible solutions: 1) a young burst or ongoing constant SF with non- 

negligible extinction or 2) a composite young -I- "old" stellar population. In the 
first case one obtains a best fit E{B — T^) ~ 0.25, or Ay ^ 0.5-1.8 at a 1 cr level. 
In consequence we obtain SFR ^ 11-41 M© yr"-"^, higher than earlier estimates, 
and we estimate a fairly high total luminosity [L ^ (1 — 4) x IO^^Lq) for this 
galaxy, in the range of luminous infrared galaxies. A Lya transmission of 23- 
90 % is estimated from our best fit models. Other properties (age, SF history) 
remain largely unconstrained. In case of composite stellar populations the SFR, 
mass, and luminosity estimate is lower. The two scenarios may be distinguishable 
with IRAC/Spitzer observations at 3.6 and 4.5 /im. 

Key words: Galaxies: high-redshift - Galaxies: evolution - Galaxies: starburst - 
Cosmology: early Universe - Infrared: galaxies 



1 INTRODUCTION 

Little is known about the stellar properties, extinction, and 
the expected intrinsic Lya emission of distant, high redshift 
galaxies. Indeed, although it has in the recent past become 
possible through various techniques to detect already size- 
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able numbers of galaxies at z > 5 (see e.g. the reviews of 
Taniguchi et al. 2003 and Spinrad 2004) the information 
available on these objects remains generally scant. For ex- 
ample, in many cases the galaxies are just detected in two 
photometric bands and Lya line emission, when present, 
serves to determine the spectroscopic redshift (e.g. Bremer 
et al. 2004, Dickinson et al. 2004, Bunker et al. 2004). Then 
the photometry is basically used to estimate the star forma- 
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tion rate (SFR) assuming standard conversion factors be- 
tween the UV restframe light and the SFR, and nothing is 
known about the extinction, and the properties of the stel- 
lar population (such as age, detailed star formation histories 
etc.) 

At higher redshift [z > 6) even less information is gen- 
erally available (but see a recent study of Eyles et al. 2005 
on two z 6 galaxies observed with HST and Spitzer). 
Many objects are found by Lya emission, but remain weak 
or sometimes even undetected in the continuum (e.g. Rhoads 
& Malhotra 2001, Kodaira et al. 2003, Cuby et al. 2003, 
Ajiki et al. 2003, Taniguchi et al. 2004). In these cases the 
Lya luminosity can be determined and used to estimate a 
SFR using again standard conversion factors. Also the Lya 
equivalent width is estimated, providing some possible clue 
on the nature of these sources. However, this has lead to 
puzzling results e.g. for the sources from the LALA survey 
which seem to show unusually large Lya equivalent widths 
that are difficult to understand without invoking exceptional 
conditions (PopIII stars?; Malhotra & Rhoads 2002, Rhoads 
et al. 2003). Given the few data available for the LALA 
sources it is fair to say that the nature of these objects, their 
stellar populations, extinction etc. remain currently largely 
unknown (cf. Dawson et al. 2004). When possible, a simple 
comparison between the UV and Lya SFR is undertaken 
providing possibly information on the Lya transmission, i.e. 
the partial absorption of Lya photons on their sight line 
through the intergalactic medium (e.g. Haiman 2002, Santos 
2004) and/or on partial Lya "destruction" processes close 
to the source (e.g. due to dust or ISM geometry; Chariot 
& Fall 1993, Valls-Gabaud 1993, Tenorio-Tagle et al. 1999, 
Mas-Hesse et al. 2003). 

Notable exceptions oi z > 4 samples for which some 
estimate of extinction is available from multi-band photom- 
etry include work on the Subaru Deep Survey (Ouchi et al. 
2004) and GOODS' data (e.g. Papovich et al. 2004). Lehnert 
& Bremer (2004) also discuss some preliminary information 
on little extinction in their z > 5 sources. Interestingly, in 
their study of a 2 = 5.34 galaxy discovered by Dey et al. 
(1998), Armus et al. (1998) find indications for significant 
reddening {Av > 0.5 mag) from analysis of the observed 
SED and from the presence of Lya emission. 

In a similar manner we will here present a consistent 
study of the stellar population properties, extinction, and 
Lya emission for two galaxies at redshift 2 > 6. For this aim 
we use two distant (2 > 6) gravitationally lensed galaxies 
for which multi-band photometry is available (detection in 
at least 3-4 bands) . Through a quantitative analysis of their 
SED, using a vast library of empirical and theoretical tem- 
plate spectra, we aim to constrain properties of the stellar 
populations, such as age and star formation (hereafter SF) 
history (burst or constant SF?) and their extinction. Fur- 
thermore by comparing the Lya emission expected from the 
stellar population constraint with the observed Lya flux we 
estimate consistently the Lya "transmission" for the indi- 
vidual sources. 

The Lya transmission and SF properties derived here 
can in principle be used to infer the ionisation fraction of 
hydrogen in the IGM at a given redshift (cf. Haiman 2002, 
Santos 2004), a key quantity of interest for the study of the 
reionisation history of the Universe (cf. review from Barkana 
& Loeb 2001). Obviously the present "exploratory" work 



will have to be extended to larger galaxy samples, and so- 
phisticated tools will probably be needed to interpret such 
results in terms of IGM properties (cf. Gnedin & Prada 
2004). However, this approach should be complementary to 
other methods probing the reionisation history by measuring 
the Gunn- Peterson optical depth observed in quasar spectra 
as a function of redshift (e.g. Becker et al. 2001, Fan et al. 
2003), or by comparing Lya luminosity functions at different 
redshifts (e.g. Malhotra & Rhoads 2004). 

The remainder of the paper is structured as follows. In 
Sect.|5|we summarise the adopted observational constraints 
from the literature. Our modeling technique is described in 
Sect. The detailed results for each galaxy are presented 
in Sects. |3 and 2] Our main conclusion are summarised in 
Sect. El 



2 OBSERVATIONAL CONSTRAINTS 

The two galaxies studied here are: 1) The probable 2^7 
galaxy recently discovered by Kneib et al. (2004, hereafter 
KESR) , which presently lacks of a spectroscopic redshift but 
for which rather accurate multi-band HST observations are 
available, allowing us in particular also to derive a fairly 
reliable photometric redshift. 2) the 2 — 6.56 Lya emitter 
HCM 6A behind the lensing cluster Abell 370. We now sum- 
marise the observational data, taken from the literature. The 
adopted redshift and gravitational magnification factors are 
listed in Table [T] 

Before proceeding let us mention for clarity that these 
two objects are generally considered to be star forming 
galaxies (starbursts), not AGN (narrow line - type II - or 
others), as no contradicting information is available so far. 
However, one must bear in mind that some of the interpre- 
tations presented below (and in the literature) may need to 
be revised, should this assumption be incorrect. 

Triple arc in Abell 2218: The observational data for 
this object, named Abell 2218 KESR hereafter, is taken from 
Kneib et al. (2004, hereafter KESR) and from Egami et al. 
(2005). The photometry from KESR includes observations 
with HST (WFPC2, ACS, NICMOS) in Veoew (undetected), 
I814W, Z850LP, and H16OW, and with NIRC/Keck in J. Sub- 
sequently, additional photometry was obtained with NIC- 
MOS/HST in the J band (FUOW), and with IRAC/Spitzer 
at 3.6 and 4.5 /im (see Egami et al.) For our computations 
(see below) we have used the appropriate filter transmis- 
sion curves. In particular, updated transmission curves were 
used for the ACS and NICMOS filters (M. Sirianni 2003, 
private communication; Sirianni et al. 2004; R. Thompson 
2003, private communication). 

Few brief comments concerning the photometry are 
needed here. First, KESR present photometry for two mul- 
tiple images (a and b). Apparently sources a and b differ in 
the Z850LP flux (with quoted errors of ± 0.05 mag) by 3.2 a, 
whereas the fluxes in the other filters agree well within 1 a. 
Differential lensing across the images together with sampling 
effects could be responsible for this small discrepancy. As we 
are interested in a global representative SED for this source, 
we have chosen to use the averaged photometric SED, the 
magnification factors being the same for the two images. 
Finally, we have also noted some apparent discrepancies be- 
tween the measurements reported in KESR and Egami et 
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al., the most important one being the Hieow flux, which is 
~ 15-20 % (3-4 (j) higher in the latter publication. These 
differences are mostly due to the use of different apertures 
on different repixeled/rescaled images (J. Richard, 2004, pri- 
vate communication). Again, this illustrates the difficulty in 
deriving reliable colors for extended arcs. To account for 
these small discrepancies and for the possible error under- 
estimate we therefore adopt a minimum photometric error 
of 0.15 mag in all filters. It is worth noting that photomet- 
ric errors translate into absolute flux calibration errors for 
fitting purposes. As we will see below, adopting the latter 
minimum errorbars significantly improves the SED fits. 

In addition to the photometry, the non-detection of the 
source with Keck LRIS spectroscopy provides an upper limit 
on the continuum flux between 9000 and 9300 A (KESR). 
This upper limit will be used as an additional constraint 
in our SED modeling. KESR also indicate a possible drop 
of the continuum below ~ 9800 A from their Keck II NIR- 
SPEC spectrum. For various reasons the reality of this spec- 
tral break is questionable. First a true neutral hydrogen 
break ("Lya break") at ~ 9800 A, far in the red wing of the 
Z850LP fllter, seems incompatible with the relatively strong 
flux measured in this fllter. Furthermore, test computations 
show that such a break is difficult if not impossible to recon- 
cile with our spectral modeling. In any case the signiflcance 
of this flnding appears questionable as the detected contin- 
uum is extremely faint and noisy. The reality of this spectral 
feature is now also questioned by Egami et al. (2005). For 
these reasons this information is discarded from our spectral 
fitting. 

In practice we have retained the following two variants 
to describe the observed SED of this source: SEDl ) The 
average fluxes (I8i4w, zssolp, Hhow, Hieow) of images a 
and b from KESR plus the IRAC/Spitzer data of image 
b from Egami et al. SED2) All fluxes from image b from 
Egami et al. These are treated as SEDs from two different 
objects. Furthermore, for each of these "objects" we have 
computed two cases in our SED fitting described below: t) 
The observed SEDs in I8i4w, zssolp, Hhow, Hieow, 3.6, and 
4.5 lira, ii) Same as (i) plus the flux limits from the Veoew 
and Keck LRIS non-detections. 

No emission line has so far been detected for Abell 2218 
KESR. Its spectroscopic redshift remains therefore presently 
unknown but the well-constrained mass model for the cluster 
strongly suggests a redshift z > 6 for this source (KESR, 
Egami et al. 2005) . The magnification factors of both images 
a and b is = 25 ± 3, according to KESR. 

Abell 370 HCM6A: The observational data of this 
z — 6.56 galaxy is taken from Hu et al. (2002). The pho- 
tometry includes VRIZJHK' from Keck I and II (LRIS 
and Echellette Spectrograph and Imager) and from Sub- 
aru (CISCO/OHS). The gravitational magnification of the 
source is = 4.5 according to Hu et al. (2002). 

Photometric fluxes and errors were adopted from their 
Fig. 3. Where possible the appropriate filter transmission 
curves were used. The "Z" band filter transmission is some- 
what uncertain, as these observations were undertaken using 
an RG850 fllter, which together with the LRIS optics and 
the CCD response, yields a transmission similar to a Z band 
filter (Hu et al. 1999). Our approximate filter curve shows a 
blueward shift of Acff by ^ 200 A compared to the informa- 
tion given by Hu et al. (1999). However, since the redshift 



-4 



PopIII: pure stellar 



SFR = const 




log(age) 



Figure 1. Temporal evolution of the UV slope f3 measured be- 
tween 1300 and 1800 A from synthesis models of different metal- 
licities and for instantaneous bursts (solid lines) and constant SF 
(long dashed lines). Black lines show solar metallicity models, 
red lines metallicities between Z = 10~^ and zero (PopIII), blue 
lines intermediate cases of Z = 0.004 and 0.0004. The dotted lines 
show 13 if nebular continuous emission is neglected, i.e. assuming 
pure stellar emission. Note especially the strong degeneracies of 
l3 in age and metallicity for bursts, the insensitivity of /3 on Z for 
constant SF, and the rather red slope for young very metal-poor 
bursts. Further discussions in the text 



is known for this source and since we adjust the observed 
flux (not magnitude) in this band, this should not affect our 
conclusions. 



3 SED MODELING 

3.1 Main restframe UV-optical SED features of 
high-z galaxies and their "information 
content" 

Before proceeding to the fits of the individual SEDs a brief 
comment on the available SED features seems appropriate. 

For obvious reasons the available SED (basically from 
broad-band photometry) of high-z {z > 6) galaxies is pri- 
marily limited to the rest-frame UV (when observed from 
the ground) or optical spectrum (when available e.g. with 
Spitzer and future satellite missions). The main information 
"encoded" in this SED is therefore: 1 ) the neutral HI break 
shortward of Lya (hereafter the "Lya" break) due to the 
strong or complete Gunn-Peterson trough, 2) the slope of 
the UV spectrum, and 3) possibly a 4000 A break (hereafter 
denoted Balmer break), if present and covered by the ob- 
servations. In addition the presence of the Lya line, mostly 
used to determine spectroscopically the redshift, provides 
clear evidence for ongoing massive star formation (hereafter 
SF). 

The position of the Lya break depends essentially on 
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redshift. The UV slope depends on the intrinsic spectrum - 
in turn depending mostly on age and SF history - and on 
the extinction, i.e. the extinction law and the amount of red- 
dening. The Balmer break becomes visible (in absorption) 
in the continuum of stellar populations after > 10-30 Myr. 
Lya emission, if due to stellar photoionisation and not AGN 
activity, indicates the presence of young (< 10 Myr) massive 
ionizing stars. 

Concerning the UV slope, it is useful to recall that this 
quantity ^ does not lend itself to determine the metallicity 
of a star forming galaxy from a theoretical point of view and 
in terms of individual objects. The reasons are that intrinsi- 
cally the UV slope shows only small variations with metal- 
licity, and that the slope depends strongly on the exact SF 
history (see e.g. Leitherer & Ifeckman 1995, Meurer et al. 
1995). This is illustrated in Fig. where (3 (measured over 
the interval 1300-1800 A) is plotted as a function of age for 
populations of metallicities between solar and zero (PopIII) 
and for the limiting cases of bursts and SFR=const. Further- 
more, as pointed out in Schaerer (2002, 2003) and also shown 
in this figure, for very low metallicities {Z < 1/50^©) neb- 
ular continuous emission becomes dominant even down to 
UV wavelengths (longward of Lya), such that the observed 
integrated (stellar -|-nebular) spectrum has even a flatter UV 
slope than high metallicity starbursts. In other words, even 
for bursts, there are strong intrinsic degeneracies of /3 be- 
tween age and metallicity, to which the additional effect of 
reddening must be added, including the uncertainties on the 
a priori unknown extinction law. It is therefore evident that 
on an individual object basis there are in general degen- 
eracies between age, metallicity, SF history, and extinction. 
However, this does not preclude the possible existence of 
statistical correlations between quantities such as e.g. (3 and 
metallicity in large samples of galaxies, as known to hold e.g. 
for local UV selected starbursts (cf. Heckman et al. 1998). 
Also, as we will see below, there are cases where the UV 
slope and the mere fact of the presence of an emission line 
allow us nevertheless to lift some degeneracies and therefore 
to determine interesting constraints on the stellar popula- 
tion and on extinction. 

The behaviour of the 4000 A break and its use as an 
age indicator has extensively been discussed in the literature 
(e.g. Bruzual 1983, and recently Kauffmann et al. 2003). 
As in simple stellar populations its amplitude is basically 
a monotonically increasing function of age, an estimate of 
the break, e.g. obtained from 3.8-4.5 /im photometry with 
IRAC/Spitzer and JHK photometry in 2; > 6 galaxies, pro- 
vides information on the age of the light emitting stellar 
population. Since the exact SF history cannot be determined 
in these cases (in contrast to studies at low z, cf. Kauffmann 
et al.) the amplitude of the break provides a range of ages, 
the minimum age being given by instantaneous bursts, the 
maximum age from models with constant SF. For obvious 
reasons this maximum "luminosity weighted" age derived 
from a measure of the Balmer break is also unaffected by 
considerations of possible multiple stellar populations. The 



^ Various definitions of the UV slope exist. The most commonly 
used ones, generally denoted /3, are defined as the power-law index 
of the SED in F\ versus A over a certain wavelength interval. 



same is not true for the minimum age, which is, however, of 
less cosmological interest. 

3.2 Spectral fitting 

For the spectral fitting we use a slightly adapted version 
of the photometric redshift code Hyperz of Bolzonella et al. 
(2000). Hyperz does standard SED fitting using a number 
of modeling parameters. The free parameters for the SED 
modeling are: 

1) the spectral template, 

2) extinction and the reddening law, 

3) a parameter /Lyf describing possible deviations from the 
average Lyman forest attenuation from Madau (1995). 

For Abell 2218 KESR the source redshift is also a free pa- 
rameter. 

For the spectral templates we use a large compilation of 
empirical and theoretical SED, including starbursts, QSO, 
and galaxies of all Hubble types, and covering various star 
formation histories (bursts, exponentially decreasing, con- 
stant SF) and various metallicities. For most applications 
we group the templates in the following way: 

• Starbursts and QSOs (hereafter SB+QSO): this 
group includes the starburst templates with E{B — V) from 
< 0.1 to 0.7 from the Calzetti et al. (1994) and Kinney et al. 
(1996) atlas, the HST QSO template of Zheng et al. (1997), 
as well as UV-optical spectrum of the metal-poor galaxy SBS 
0335-052 with numerous strong optical emission lines (and 
an extinction of E{B - F) - 0.09, Izotov & Thuan 1998) 
kindly communicated to us by Yuri Izotov (2002, private 
communication) 

• BCCWW+: Bruzual & Chariot (1998, private com- 
munication; cf. Bruzual & Chariot 1993) evolving synthe- 
sis models assuming bursts, constant star formation, and 
exponentially decaying star formation histories reproducing 
present day spectra of galaxies of various types (E, SO, Sa, 
Sb, Sc, Sd, and Im) plus the empirical E, Sbc, Scd, and Im 
templates from Coleman et al. (1980), as included in the 
public Hyperz version. 

• S03+: Theoretical templates of starburst galaxies from 
Schaerer (2003) covering metallicities of Z = 0.02 (solar), 
0.008, 0.004, 0.001, 1/50 Zq, Z = 10"^ 10"^ and zero 
metallicity (PopIII). For low metallicities {Z ^ 10~^) these 
templates have been computed for 3 different assumptions 
on the IMF. The spectral library includes burst models and 
models with a constant star formation rate (SFR). For more 
details see Schaerer (2003). For the present work these com- 
putations were extended to cover ages of up to 1 Gyr. These 
SEDs are available on request from the first author and on 
the Web^. 

The standard extinction law adopted here is the one 
from Calzetti et al. (2000) determined empirically from 
nearby starbursts. We also explore the possible implications 
of other laws, such as the Galactic law of Seaton (1979) in- 
cluding the 2200 A bump, and the SMC law from Prevot et 
al. (1984) and Bouchet et al. (1985) showing no UV bump. 
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but a steeper increase of the extinction in the UV compared 
to Calzetti et al.. 

For the Lyman forest attenuation, Hyperz follows 
Madau (1995). However, we allow for possible deviations 
from the mean attenuation by varying the Lyman forest op- 
tical depths T^g^ by a multiplicative factor taking the values 
of (/Lyf, 1-, and 1//Lyf). Typically we adopted fhyi = 2 or 
3. Here r^g^ stands for the optical depths corresponding to 
the absorption between Lya and Ly/3, and between Ly/3 and 
the Lyman limit respectively. 

The following other minor changes have been made in 
our version (1.3ds) of Hyperz. The calculation of the syn- 
thetic photometry deals correctly with templates including 
strong spectral lines (emission or absorption). Furthermore 
we make sure to use the proper filter transmission curves 
usually given in photon units. Earlier versions of Hyperz and 
other codes (e.g. evolutionary synthesis codes) assume some- 
times (for "historical" reasons) that transmission curves be 
given in flux units. In case of wide filters, e.g. such as some 
ACS/HST filters, this may lead to small differences. Other 
modifications concern essentially features related to the user 
interface (additional outputs etc.). 

For given choices of the above parameters, theHyperz 
code performs a standard minimisation fit to the observed 
SEDs and determines, for each point in the parameter space, 
the corresponding value. Using these values, it is pos- 
sible to quantify the probabilities for the main free param- 
eters, namely extinction, age of the spectral template, SF 
history, etc. When the SED fitting is based on theoretical 
templates, the SFR value is easily obtained and allows us to 
compare the expected values for the Lya flux to the actual 
ones. 

To convert the observed/adjusted quantities to abso- 
lute values we adopt the following cosmological parameters: 
ilm = 0.3, = 0.7, and Ho = 70 km s~^ Mpc"\ 



Abell 2218 KESR - all filters with detections 
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Figure 2. Photometric redshift probability distributions P{z) 
of Abell 2218 KESR using three spectral template groups. Solid 
line: BCCWW+ template group, dotted: SB-l-QSO, long dashed: 
S03+. The three upper blue curves stand for the average SED 
of a and b (SEDl), the lower red ones for the SED of object b 
(SED2) from Egami et al. (2005). In all cases a minimum photo- 
metric error of 0.15 mag was adopted. The P{z) shown here has 
been computed based on all filters in which the object is detected 
(I814W to 4.5 Atm). 



4 RESULTS FOR ABELL 2218 KESR 

4.1 Photometric redshift estimate 

As a spectroscopic redshift has not been obtained (yet) for 
this galaxy we here examine its photometric redshift esti- 
mate. In Fig.|21we show the photometric redshift probability 
distributions P(z) for the two SEDs (SEDl, SED2) of Abell 
2218 KESR described above using the three spectral tem- 
plate groups and adopting a minimum photometric error of 
0.15 mag. For each redshift, P{z) quantifies the quality of 
the best fit model obtained varying all other parameters (i.e. 
extinction, /Lyf, spectral template among template group). 
Given the excellent HST (WFPC2, ACS and NICMOS) pho- 
tometry, P{z) is quite well defined: the photometric redshift 
ranges typically between Zphot ~ 5.5 and 7.3. Outside of the 
plotted redshift range P{z) is essentially zero. If we assume 
the (smaller) quoted formal photometric errors (but note 
the discrepancies discussed in Sect. |5J P{z) becomes more 
peaked, i.e. the photometric redshift better defined. This is 
driven by the error on the zgsoLP flux, which determines the 
red side of the "Lya" break. However, the resulting best 
fit value 2phot does not change much. Furthermore, the fit 
quality is considerably decreased. This demonstrates the in- 
terest of such high accuracy measurements and the need for 
reliable error estimates. 



The predicted redshift distribution is found to be quite 
insensitive to the exact template (as shown in Fig. to 
the exact value of /Lyf, and to the adopted extinction law 
(variations of the latter two are not shown). However, we 
note that for this object the fits (and P{z)) are improved 
when allowing for deviations from the average Madau (1995) 
attenuation law. The curves shown here have been computed 
for /Lyf = 2. 

More important in determining P{z) is the exact SED. 
As seen from Fig.Hthe use of SEDl or SED 2 lead to some- 
what different P{z) distributions. SED2 (cf. Sect. |2l yields 
a somewhat larger redshifts, albeit with a somewhat yre- 
duced fit quality. These differences illustrate how uncertain- 
ties and difficulties in the photometric measurements of such 
faint sources, whose origin are briefly discussed in Sect. ^ 
propagate to the photometric redshift estimate. 

All our best-fit solutions have redshift Zphot '~ 6.25- 
6.63, lower than the redshift range estimated by KESR, but 
compatible with the more recent quantitative analysis of 
Egami et al. (2005). Given the various free parameters, un- 
certainties on the intrinsic SED, etc. we conclude that the 
redshift of AbeU 2218 KESR is likely a - 6.0 - 7.2. taking 
into account both our photometric determination and the 
lensing considerations of KESR. 
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Figure 3. Best fits SEDs to the observations of Abell 2218 (SED2 
from Egami et al. 2005, including the flux limit from the non- 
detections in Vgoew and at 9000-9300 A from spectroscopy; cf. 
Sect.|3. The red crosses indicate the corresponding model broad 
band fluxes. The solid line shows the best fit for a template from 
the S03-I- group, and dotted from the SB+QSO group. The red- 
shift for these solutions are z ~ 6.63 and 6.54 respectively. See 
text for more information 

4.2 SED fits and inferences on the stellar 
population and on Lya. 

A large number of models have been computed using the 
different variants of the SEDs describing this object (SEDl- 
2), the different filter combinations (non-detections 4- spec- 
troscopic constraint) discussed in Sect. |51 and varying the 
various model parameters. We first discuss briefly the main 
salient results with the help of some illustrations. A more 
general discussion of the results and their dependence on 
various assumptions follows. 

4-2.1 Age, star formation history, and extinction 

Figure |3 shows the best fit models to the SED2 including 
the upper limits from the I8i4w and LRIS spectroscopy for 
the S03-I- and SB+QSO template groups. The best fit red- 
shifts are Zphot = 6.63 and 6.54 respectively. These fits show 
in particular that the spectroscopic constraint can be ac- 
commodated simultaneously with the observed zgsoLP flux; 
the resulting fits are within the 1 a errors in all bands. The 
best fit from the S03-(- group corresponds to a burst with an 
age of 15 Myr at solar metallicity and no extinction (solid 
line). Similarly good fits are also obtained for lower metallic- 
ity. The best fit with empirical starburst and QSO templates 
(SB-I-QSO) is obtained with the spectrum of the metal-poor 
H II galaxy SBS 0335-052 (dotted fine in Fig. EJ. In this 
case the apparent Balmer break observed between the NIC- 
MOS/HST and IRAC/Spitzer domain is simply explained 
by the presence of strong emission lines in the 3.6 and 4.5 



Figure 4. Best fits SEDs to the observations of Abell 2218 (Large 
symbols: SED2 from Egami et al. 2005; small symbols: SEDl; 
cf. Sect. Only true detections are taken into account. The 
red crosses indicate the corresponding model broad band fluxes. 
The solid (dotted) line shows the best fit for a template with 
constant SFR from the S03+ group to the SEDl (SED2). The 
best "maximum age" fits correspond to ages of 500 and 400 Myr, 
no extinction, and redshifts Zphot ~ 6.40 and 6.57 respectively for 
SEDl and SED2. See text for more information 

/xm filters ^ and some additional extinction to reduce the 
restframe UV fiux. The extinction needed is Ay = 0.6 for 
the Calzetti et al. law, or Ay ~ 0.2 for the Prevot et al. 
extinction law. In terms of age the restframe UV to opti- 
cal spectrum (continuum and lines) of SBS 0335-052 corre- 
sponds to a young population of ~ 3-5 Myr according to the 
analysis of Papaderos et al. (1998) and Vanzi et al. (2000). 
Of course, the presence of an older population in addition 
to the starburst cannot be excluded on the present grounds. 
In short, the observed SED of Abell 2218 KESR can be 
explained by a young population withour or with emission 
lines. Spectroscopy in the 3-4 /im range would be needed to 
distinguish the latter solution from others. 

Alternatively, good SED fits are also obtained with rela- 
tively "old" populations. The oldest ages are obtained when 
invoking the longest SF timescale, i.e. constant SF. In this 
case the UV restframe fiux remains high (due to the contin- 
uous formation of massive stars) and older ages need to be 
attained to build up a sufficient population of evolved stars 
with strong Balmer breaks, in order to reproduce the ob- 
served break. This case is shown in Fig.|l]with fits with the 



^ The main lines are between H7, H/3 and [O in] AA4959,5007 
in the 3.6 fira filter and He I A5876 in the 4.5 fim filter. E.g. for 
the emission lines between H7 and [O in] AA4959,5007 the total 
observed equivalent width (boosted by the (1 -|- 2) factor) is ~ 
9130 A, as estimated from the data of Izotov & Thuan (1998), 
compared to a filter surface of 6600 A. 
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Ch\2 contours - Abell 2218 KESR (z = 6.55, S03 + , SFR = const) 
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age [Gyr] 

Figure 5. contour plot in extinction - age for solutions fitting 
the observed SED2 at redshift z = 6.55 with templates from 
the S03+ group assuming constant SF. The best fit solution is 
indicated by the black dot. Equidistant levels with a spacing of 
1 are shown. The (ID) 68, 90, and 99 % confidence regions (Ax^ = 
1., 2.71, 6.63) are delimited by the thick black lines in long dashed, 
dotted and solid respectively. Corresponding formation redshifts 
ZfoTui (assuming the cosmological parameters given in Sect. 13.21 
are indicated by the arrows. Discussion in text 



S03+ templates to the observed SEDl and SED2. The best 
fits solutions obtained here correspond to ages of 500 and 
400 Myr, no extinction, and redshifts Zphot ~ 6.40 and 6.57 
respectively. Similar, but somewhat older ages are obtained 
for metallicities Z < 0.008, below the ones shown here. 

A quantitative examination of the "maximum age" al- 
lowed by the observations (here SED2) is presented in Fig. 
13 showing contours in the extinction-age plane for a 
given set of spectral templates (S03-I- group with constant 
SFR), the Calzetti et al. extinction law, and a fixed redshift 
of z = 6.55. For these conditions the best fit is corresponds 
to 400 Myr, zero extinction, and Zphot ~ 6.57 (see Fig. IIJ. 
This Figure shows that a maximum age of ~ (250-650) Myr 
(1 <j interval) is obtained, in good agreement with the mod- 
eling of Egami et al. (2005). If true, this would correspond 
to a formation redshift of Zform ~ 8.7-20 for our adopted 
cosmological parameters. As clear from Fig. |K| even with 
constant SF models, younger populations with some extinc- 
tion can also fit, although less well, the present observations. 
However, solutions with low or zero extinction are generally 
preferred. 

When varying the star formation history between these 
extreme cases (burst or SFR=const), i.e. considering e.g. ex- 
ponentially declining SF histories, any intermediate age can 
be found for obvious reasons. Such cases are e.g. obtained 
when fitting templates from the Bruzual & Chariot models 
(not shown here) with exponentially decreasing SF histo- 
ries and can be found in Egami et al. As discussed in Sect. 



13.11 considering multiple stellar populations (cf. Eyles et al. 
2005) does not alter the above estimate of the maximum 
age determined from constant SFR models. In any case, the 
data available here does not allow us to constrain the SF 
history and age further. 

4-2.2 Lya emission 

The observations obtained so far have not revealed any emis- 
sion line from this object (KESR). In particular Lya emis- 
sion is lacking, which could be puzzling for a source with 
intense star formation. As we have just seen a variety of 
star formation histories and ages are possible for Abell 2218 
KESR. Therefore one may or may not expect intrinsic Lya 
emission. 

A simple explanation for the apparent absence of Lya 
emission could be to invoke an advanced age (in the post 
starburst phase). However, even with a young age it is not 
necessary that Lya emission be observed. E.g. the spectrum 
of the metal-poor H ll galaxy SBS 03335-052 which provides 
an excellent fit to the observed SED and shows strong emis- 
sion lines (cf. Fig. ^ does not show Lya emission (Thuan 
et al. 1997, Kunth et al. 2003). Alternatively, if intrinsically 
present, the Lya non-detection could be due a variety of fac- 
tors: a redshift z < 6.4 placing Lya below the spectral range 
discussed in detail by KESR*, a flux below their strongly 
varying detection threshold, or other factors depressing the 
Lya emission within the host galaxy (dust, ISM-f-HI geom- 
etry) and in the intervening IGM. 

In conclusion, from the available data the apparent lack 
of Lya emission from this source is not puzzling. However, 
it is not completely excluded that the galaxy truely shows 
Lya emission, which has so far eluded detection. 

4-2.3 General comments on fits and discussion 

After these main findings we shall now quickly mention more 
"technical" results about the infiuence of various fit param- 
eters. 

Quite generally, the results on the age, SF history, ex- 
tinction etc. depend little on the different variants of the 
observed SEDs (SEDl-2), on the inclusion or not of the non- 
detections in the fits, and on the use of the published formal 
errors or a minimum error of 0.15 mag (cf. Sect.|2l. The 
results discussed above are therefore quite robust with re- 
spect to these assumptions. Small differences in the best fit 
values can, however, be obtained. E.g. the best fit photomet- 
ric redshift can vary by up to < 0.2 depending on adopting 
SEDl or SED2. In all cases we note that SEDl allows bet- 
ter fits (smaller x^) than SED2. Adopting (Tmin ~ 0.15 mag 
also significantly increases the fit quality. Finally, consider- 
ing variations around the mean Lyman forest attenuation 
improves the fits (especially as the HST photometry deter- 
mining the "Lya break" is quite accurate). In practice all 
best fits are found with an increased Lyman-forest opacity 

(/Lyf =2). 

* This is probably excluded as, according to J. -P. Kneib (2005, 
private communication), no emission line was found in the blue 
part of the spectrum taken with LRIS. 
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To summarise, given the absence of a spectroscopic red- 
shift, a fair number of good fits is found to the observations 
of Abell 2218 KESR when considering all the free parame- 
ters. The main conclusions from these "best fits" are: 

1 ) Generally the determined extinction is negligible or zero 
quite independently of the adopted extinction law. The best 
fit with the empirical starburst spectrum of SBS 0335-052 
represents an exception to this case, requiring an additional 
Av ~ 0.2-0.6 mag, depending on the adopted extinction 
law. 

2) Although generally burst models fit somewhat better 
than those with constant star formation among the theo- 
retical templates (BC, S03-|~), the data does not strongly 
constrain the star formation history. 

3) Typical ages between ~ 15 and 400 Myr are obtained. 
A reasonable l-u upper bound on the age of ~ 650 Myr can 
be obtained assuming constant star formation. However, the 
data can also be well fit with a very young (~ 3-5 Myr) 
stellar population with strong emission lines (using e.g. the 
spectrum of the metal-poor galaxy SBS 0335-052). In this 
case the apparent Balmer break observed between the HST 
and Spitzer broad-band photometry is simply due to the 
presence of strong emission lines affecting the red 3.6 and 
4.5 ^m filters. 

4) Given degeneracies of the restframe UV spectra be- 
tween age and metallicity (cf. above) no clear indication on 
the galaxian metallicity can be derived, in contrast to the 
claim of KESR. Good fits to the available data can even be 
found with solar metallicity starburst templates. 

5) Depending on the star formation history and age one 
may or may not expect intrinsic Lya emission, i.e. an im- 
portant H II region around the object. The apparent absence 
of observed Lya emission does therefore not provide much 
insight. 

A more complete error analysis beyond the level pre- 
sented here is difficult to achieve for a variety of reasons 
and clearly beyond the scope of this publication. 



4-2. i SFR, stellar mass and luminosity 

The theoretical templates can also be used to estimate the 
stellar mass involved in the starburst or the star formation 
rate when constant star formation is assumed. For this aim 
we use all the best fits to the three SEDs (SEDl-3) with 
the S03-I- templates, we assume a typical redshift of z = 
6.6, and the magnification /i = 25 determined by KESR. 
For the adopted cosmology the distance luminosity is then 
dL =64457.8 Mpc. 

When constant SF is assumed one obtains the follow- 
ing star formation rate; SFR ~ (0.9 — 1.1) Mq yr~^ (for 
a Salpeter IMF from 1 to 100 M©). For the best fit ages 
of ^ 400-570 Myr the total mass of stars formed would 
then correspond to ~ (3.6 — 6.3) x 10* Mq. The mass es- 
timated from best fit burst models (of ages ~ 6-20 Myr) 
is shghtly smaller, M^, ~ (0.3 - 1) x 10* Mq. If we assume 
a Salpeter IMF with Miow = 0.1 Mq the mass and SFR 
estimates would be higher by a factor 2.55, and in good 
agreement with the values derived by KESR and Egami et 
al. (2005). In all the above cases the total luminosity (un- 
lensed) is typically Lboi ~ 2 x 10^" Lq. 



5 RESULTS FOR ABELL 370 HCM 6A 

5.1 SED fits and inferences on the stellar 
population 

Overall the pubhshed SED of HCM 6A (see Fig. [SJ is "red- 
dish" , showing an increase of the fiux from Z to H and even 
to K'^. From this simple fact and the above explanations 
it is already clear qualitatively that one is driven towards 
solutions with a) "advanced" age and little extinction or b) 
constant or young star formation plus extinction. However, 
a) can be excluded as no Lya emission would be expected 
in this case. 

Quantitatively, the best solutions obtained for each 
spectral template group is shown in the left panel of Fig- 
ure |S| Indeed, the solutions shown correspond to bursts of 
ages ~ 50-130 Myr (BCCWW+, S03 templates) and little or 
no extinction. However, as just mentioned, solutions lacking 
young (< 10 Myr) massive stars can be excluded since Lya 
emission is observed. The best fit empirical SB-I-QSO tem- 
plate shown corresponds to the spectrum of the H ll galaxy 
SBS 0335-052 with an additional extinction of Ay = 1. To 
reconcile the observed SED with Lya, a young population, 
e.g. such as SBS 0335-052, or constant SF is required. In 
any of these cases fitting the "reddish" SED requires a non 
negligible amount of reddening. 

To illustrate the typical range of possible results we 
show in Fig. Q contour maps and corresponding confi- 
dence intervals for solar metallicity models (S03-I- template 
group) and reddened with the Calzetti law. The left panel 
(burst models) illustrates in particular the need for progres- 
sively higher extinctions the younger the bursts. All ages > 
10 Myr are excluded for the absence of Lya emission. From 
the constant SF models (right panel) we see that for a given 
age Av is typically ~ 0.5-1.8 mag at the 68 % confidence 
level. For obvious reasons, no constraint can be set on the 
age since the onset of (constant) SF. 

Hence, from the photometry of HCM 6A and from the 
presence of Lya we are led to conclude that this object must 
suffer from reddening with typical values of Ay ~ 1. for 
a Calzetti attenuation law. A somewhat smaller extinction 
{Av ~ 0.4) can be obtained if the steeper SMC extinction 
law of Prevot et al. (1984) is adopted. From the present 
data it is not possible to distinguish the different extinc- 
tion/attenuation laws. Also, it is not possible to draw any 
constraints on the metallicity of HCM 6A from the available 
data (cf. Sect.lOl. 

What if we are dealing with composite stellar pop- 
ulations? Indeed it is conceivable that the Lya emission 
originates from a population of young stars and the "red- 
dish" restframe UV fiux be due to another, older popula- 
tion. Assuming constant SF, no loss of Lya and standard 
SFR conversion factors, the observed Lya emission implies 
a maximum UV fiux of the order of 0.1 /iJy (and approxi- 
mately constant in F,j over the observed wavelength range: 
Areat ^ 3000 A) for an unreddened population. The bulk of 
the observed flux could then be from an older population. 
In this case the rising spectrum from the zgsoLP over the JH 
(and presumably K) bands could even be due to an unred- 

^ The significance of a change of the SED slope between ZJH 
and HK' seems weak, and difficult to understand. 




Figure 6. Best fits SEDs to the observations of Abell 370 HCM 6A. The red crosses indicate the corresponding model broad band 
fluxes. Solid lines show the best fit for a template from the BC+CWW group, dotted from SB+QSO group, and dashed from the S03+ 
group (see explanations in Sect.^J. Left: Observed spectral range. Right: Predicted SED in Spitzer/IRAC domain for best fit models. 
Dashed lines show the bursts from the BCCWW+ and S03+ template groups. The dotted line is the spectrum of SBS 0335-052 from the 
SB+QSO group with additional Ay = 1. The solid lines show best fits for constant star formation using different extinction/attenuation 
laws (Calzetti starburst law versus SMC law). The solid triangles illustrate the IRAC point-source sensitivity (f a) for low and medium 
backgrounds excluding "confusion noise". 



Chi2 contours (z = 6,56 galaxy: Hu et al. 2002) Chi2 contours (z = 6,56 galaxy: Hu et al. 2002) 
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Figure 7. contour plots showing solutions in extinction - age diagrams. The best solutions are indicated by the black dot. Equidistant 
levels with a spacing of 0.5 are shown. The 2D 68% confidence region (corresponding to Ax^ = 2.3) is delimited by the solid thick 
black line. The (ID) 68 % confidence region for Ay (Ax^ = 1.) at each given age is delimited by the dashed thick black line. Left: 
Plot for solutions using a solar metallicity burst template from the S03+ template group and the Calzetti attenuation law. Although 
providing a good fit to the photometry, the region corresponding to ages > 15 Myr, right of the dotted vertical line, is excluded as no 
emission line would be expected in this case. Right: Same as left panel for constant star formation models. The solutions indicate a 
non-negligible extinction, but no constraint on age. Discussion in text. 
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Table 1. Summary of the main adopted and estimated properties of the analysed high redshift galaxies. The adopted magnification /i 
is given in col. 2, col. 3 gives the redshift, col. 4 an indication on the most plausible star formation history (burst or constant SF), col. 
5 a plausible age of the stellar population, col. 6 an estimate of the optical extinction Ay (for the Calzetti et al. 2000 law), col. 7 the 
estimated SFR (for a Salpeter IMF from 1-100 Mq), col. 8 an estimated stellar mass (for same IMF), col. 9 the estimated bolometric 
luminosity, and col. 10 the estimated Lya transmission (fraction of intrinsic emitted flux over observed Lya flux). See Figs . 171 and ISl and 
the corresponding text for an estimate of the confldoncc levels and range of parameters. 



Object 




redshift 


SF history 


age 
[Myr] 


Av 
[mag] 


SFR 

[Moyr-i] 


stellar mass 
[Mq] 


Lhoi 
[L©] 


Lya 

transmission 


Abell 2218 
KESR 


25. 


~ 6.0-7.2 


? 


3-400 


negligible? 


~ 1 


(0.3-6) X 10* 


~ 2 X IQi" 




Abell 370 
HCM6A 


4.5 


6.56 


CSFR/ 
young burst 


? 


~ 1. 


11-41 


(1-4) xlO* 


(1-4) xlO" 


23-90 % 








composite 


young+ "old" 


negligible 


> 0.4-0.8 


~ 2 X 10^ 


~ 3 X 10^" 


>40 % 



dened population; a strongly increasing flux and probably a 
significant "Balmer" break is then expected, similar to the 
aged burst shown in the right panel of Fig. |S] This explana- 
tion should in principle be testable with Spitzer observations 
as discussed below. 

How does our possible indication for a high extinction 
fit in with other studies? At redshift z < 4 the extinction of 
Lyman break galaxies (LBGs) has been estimated by various 
authors (e.g. Sawicky & Yee 1998, Meurer et al. 1999, Adel- 
berger & Steidel 2000, Shapley et al. 2001, Ouchi et al. 2004). 
Given their mean/median values (typically < E{B — V) >~ 
0.15-0.2) and the E{B - V), our finding of a "high" extinc- 
tion is not exceptional. Furthermore, Armus et al. (1998) 
find indications for Ay > 0.5 mag from their analysis of a 
z — 5.34 galaxy. However, this implies that dust extinction 
is likely present in starburst galaxies with redshifts above 6. 
Large amounts of dust have already been observed in QSOs 
up to similar redshift (Bertoldi et al. 2003, Walter et al. 
2003). 

5.2 Properties of HCM6A: SFR, mass, Lya 
transmission 

To estimate properties such as the mass, SFR, and the in- 
trinsic Lya emission from HCM 6A we simply examine the 
predictions from the best fit models, scale them appropri- 
ately to the luminosity distance^, and correct for the gravita- 
tional magnification (here jj, = 4.5). The derived quantities 
are summarised in Table Q 

First we consider single (non-composite) stellar popu- 
lations. From the best fit constant SF models (with variable 
ages) we deduce an extinction corrected star formation rate 
of the order of SFR(UV) ~ 11 - 41 Mq yr'^ for a Salpeter 
IMF from 1 to 100 Mq. For a commonly adopted, although 
unjustified, Salpeter IMF down to 0.1 Mq this would in- 
crease by a factor 2.55. Actually this estimate is not very 
different than the one obtained from standard SFR calibra- 
tions provided the same assumptions on the IMF. Indeed 
the observed restframe UV luminosity, e.g. derived from the 
average J, H, and K' flux of Ftjv = (2.6 ± 0.7) x 10"^° erg 

^ For the adopted cosmology and z = 6.56 one has = 64005.7 
Mpc. 



s"^ cm-2 Hz"^ is Lrcstuv = iirdlFw + z)/ /i ?a 4. x 10^** 
erg s"^Hz"^ (Hu et al. 2002), translates to SFRw 
c Lrostuv lO" '*'*"^ , where c is the usual SFR conversion co- 
efficient, and ^uv the UV extinction. For the standard value 
c = 1.4 X lO"'^** from Kennicutt (1998), assuming a Salpeter 
IMF down to 0.1 Mq, and Auv 2.-3., one has SFR ~ 
35-88 Mq yr"^; for the IMF used in this work (Salpeter 
from 1-100 Mq) this becomes SFR - 14-34 Mq yr"^ This 
assumption, and the absence of extinction correction also ex- 
plains the difference with SFR estimate of Hu et al. (2002) 

7 

For continuous SF over timescales tsp longer than 10 
Myr, the total (bolometric) luminosity output is typically 
- 10^° Lq per unit SFR (in Mq yr"^) for a Salpeter IMF 
from 1-100 Mq, quite independently of metallicity. The to- 
tal luminosity associated with the observed SF is therefore 
L ~ (1 — 4) X lO^^I/Q , close to or just above the limit to possi- 
bly qualify as a luminous infrared galaxy (Lir > 10^ '^Lq; cf. 
Sanders & Mirabel 1996) if a signiflcant fraction of its bolo- 
metric flux emerges in the (restframe) IR. For tsF ~ 10 Myr 
the estimated stellar mass is « tsF x SFR ~ (1 — 4) x 10* 

Mq. 

From the data given by Hu et al. (2002), the observed 
Lya flux is F(Lya) = ^ x 3.8 x 10"^* erg s~^ cm~^, with 
the magniflcation factor /j.. The Lya luminosity per unit 
SFR from the same S03 models used above is L(Lya) = 
(2.4 - 4.4) X 10"^ erg s"^ (Mq yr"^)"^ for metallicities be- 
tween solar and 1/50 Zq. The SFR deduced from Lya would 
then be SFR(LyQ) ~ 0.4-0.8 Mq yr"^ for HCM 6A. Taking 
an extinction of Ay = 1. (for the Calzetti law) into ac- 
count implies a reddening corrected SFR(Lya) ~ 7-12 Mq 
yr"\ The ratio between SFR(Lya)/SFR(UV) presumably 
reflects the incomplete Lya transmission tLya, which can 
be estimated in various ways. The most consistent estimate 
is obtained from the comparison of the predicted Lya lu- 
minosity of each best fit model (obtained from fitting the 
broad-band SED) to the observed Lya luminosity. From the 
best fit models with Av ~ 1 and the Calzetti law we ob- 
tain ti^ya ~ 23-54 %; in the case of the fit with the Prevot 

^ Actually Hu et al. (2002) derive without further explanation 
SFR = 9 Mq yr~^ from Lrestuv = 4. X 10^* erg s~^, whereas 
the classical Kennicutt (1998) calibration would yield SFR = 5.6 
Mq yr~^ without extinction correction. 
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et al. extinction law we find a higher transmission thya '~ 
90 %. For comparison Haiman (2002) assumed a Lya trans- 
mission of 20 % from the data of Hu et al. Per definition 
this Lya "transmission" corresponds to the ratio of the ex- 
pected/intrinsic Lya emission from the starburst over the 
observed one. The physical causes for a partial (i.e. < 100 
%) transmission are of course open to various interpreta- 
tions (e.g. physical processes destroying Lya photons in the 
host galaxy, absorption in the intervening IGM etc.). In fact 
the relatively high Lya transmission estimated here could 
be somewhat surprising, given the Gunn-Peterson trough 
observations in z > 6 quasars (cf. Becker et al. 2001, Fan et 
al. 2003) and the possible presence of dust (this work). 

Consider now the case of composite stellar popula- 
tions. In this case we retain as a rough estimate in Table 
the SFR{Lya) (from the young population) as a lower 
limit, and the mass and total luminosity is derived from 
the best fit burst model of age ~ 100 Myr assuming that 
this "older" population dominates the observed continuum 
flux. Formally we then have no handle on the Lya trans- 
mission, except that it cannot be very low (say < 40 % 
« 0.1/0.26 = Fyoung/ < Fobs >) since otherwise the as- 
sociated UV flux from the young population Fyoung would 
dominate the observed continuum flux < fobs >■ 



5.3 Spitzer Observatory predictions 

It is interesting to examine the SEDs predicted by the var- 
ious models at longer wavelengths, including the rest-frame 
optical domain, which is potentially observable with the sen- 
sitive IRAC camera onboard the Spitzer Observatory and 
other future missions. In the right panel of Fig. |S|we plot 
the 3 best fits to the observed data for the BCCWW-I- and 
S03-I- template groups ( "burst" solutions with no extinction) 
and the SBS 03352-052 template (with additional Av = 1) 
showing strong optical emission lines. We see that these so- 
lutions have fluxes comparable to or above the detection 
limit of IRAC/Spitzer On the other hand the strongly 
reddened constant SF or young burst solutions do not ex- 
hibit a Balmer break and are hence expected to show fluxes 
just below the IRAC sensitivity at 3.6 /im and significantly 
lower at longer wavelengths. As Lya emission is expected 
only for the reddened SEDs the latter solutions (low 3.6-4.5 
/xm flux) are predicted to apply to HCM 6A, except if com- 
posite stellar populations are considered. Indeed, a high 3.6 
and 4.5 fim flux could be a good indication for a compos- 
ite stellar population, as discussed above. In addition it is 
important to secure higher accuracy photometry especially 
in the near-IR (JHK) to assess the accuracy of the redward 
increasing shape of the spectrum (in F^), which drives one 
towards solutions with non-negligible extinction. 



* The IRAC detection limits plotted here correspond 
to the values given by the Spitzer Science Center on 
jhttp : //ssc . spitzer . caltech. edu/irac/sens .html as 1 ct 
point-source sensitivity for low and medium backgrounds for 
frame times of 200s and described by Fazio et al. (2004). These 
values do not include "confusion noise". 



6 CONCLUSION 

Using SED fitting techniques considering a large number 
of parameters (mainly a vast library of empirical and the- 
oretical template spectra, variable extinction and extinc- 
tion laws) we have attempted to constrain the properties 
of the stellar populations and Lya emission of two strongly 
lensed galaxies with redshifts z > 6 from their observed 
SED including various ground-based observations, HST, and 
Spitzer observations. 

The following main results have been obtained for these 
objects (see Sects. |3 21 1 ^-nd summary in Table 0: 

• Triple arc in Abell 2218 discovered by Kneib et al. 
(2004, KESR). The most likely redshift of this source is 2 ~ 
6.0-7.2 taking into account both our photometric determi- 
nation and lensing considerations. 

SED fits indicate generally a low extinction {E(B — V)< 
0.05) but do not strongly constrain the SF history. Best fits 
have typical ages of ~ 3 to 400 Myr. A reasonable maxi- 
mum age of (250-650) Myr (1 g interval) can be estimated. 
However, the apparent 4000 A break observed recently from 
combination of IRAC/Spitzer and HST observations, can 
also equally well be reproduced with the template of a young 
3-5 Myr) burst where strong restframe optical emission 
lines enhance the 3.6 and 4.5 ^ra fiuxes. 

The estimated SFR is typically ~ 1 M© yr~^ for a Salpeter 
IMF from 1-100 M©, in agreement with previous estimates. 

Given the poor constraint on age and SF history, we con- 
clude that intrinsic Lya emission may or may not be present 
in this galaxy. The apparent non-detection of Lya by KESR 
can therefore even be understood without invoking Lya de- 
struction. 

• Abell 370 HCM 6A discovered by Hu et al. (2002). 
The relatively red SED and the presence of Lya emission 
indicate basically two possible solutions: 1) a young burst or 
ongoing constant SF with non-negligible extinction [Av ~ 
0.5-1.8 at a 1 (T level) or 2) a composite young -f- "old" 
stellar population. 

For the first case, best fits are obtained for constant SF 
with E{B — V)^ 0.25. In consequence previous SFR esti- 
mates for this source must likely be revised upward. If cor- 
rect, the bolometric luminosity of this galaxy is estimated 
to be 1/ ~ (1 — 4) X 10^^1/0, comparable to the luminosity of 
infrared luminous galaxies. Furthermore a Lya transmission 
of ~ 23-90 % is estimated from our best fit models. 

Alternatively the observed 0.9-2.2 SED could also be 
fit without extinction by a composite "young" and "old" 
stellar population, where the former would be responsible 
for the Lya emission and a fraction of the restframe UV flux. 
The SFR, stellar mass, and total luminosity are then lower 
than in case 1. The two scenarios may be distinguishable 
with IRAC/Spitzer observations at 3.6 and 4.5 /im. 

Given the limited observed spectral range, the present 
data does not allow to draw any firm constraints on the 
maximum age of the stellar population. 

In general it should also be noted that broad-band SED 
fits or measurements of the UV slope do not allow one to 
determine the metallicity of a star forming galaxy from a 
theoretical point of view and in terms of individual objects 
given important degeneracies (cf. Sect. 13.111 . 

The estimates of the Lya transmissions presented here 
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can in principle be used to constrain the intervening IGM 
properties, and therefore probe the reionisation of the Uni- 
verse. 

Although the results obtained here from this ex- 
ploratory study of just two lensed galaxies, the highest 
known redshift galaxies with photometric detections in at 
least 3-4 filters, cannot provide a general view on the SF 
and IGM properties at z > 6, there is good hope that the 
sample of such objects will considerably increase in the near 
future with the availability of large ground-based telescopes 
and sensitive space borne observatories such as Spitzer and 
even more so with the planned JWST. 
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